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fRACT 


A FORTRAN computer code (SKINTEMP) has been developed to calculate 
transient missile/aircraft aerodynamic heating parameters utilizing basic flight 
parameters such as altitude, Mach number, and angle of attack. The insulated skin 
temperature of a vehicle surface on either the fuselage (axisymmetric body) or 
wing (two-dimensional body) is computed from a basi'* ’.eat balance relationship 
throughout the entire spectrum (subsonic, transonic, _>i ~ .« sonic, hypersonic) of 
flight. This calculation method employs a simple Unite difference procedure which 
considers radiation, forced convection, and non-reactive chemistry. Surface 
pressure estimates are based on a modified Newtonian flow model. Eckert’s 
reference temperature method is used as the forced convection heat transfer model. 
SKINTEMP predictions are compared with a limited number of test cases. 
SKINTEMP was developed as a tool to enhance the conceptual design process of 
high speed missiles and aircraft. Recommendations are made for possible future 


development of SKINTEMP to further support the design process. 







THESIS DISCLAIMER 


The reader is cautioned that the computer program SKINTEMP 
has not been exercised for all cases of interest. While every 
effort has been made within the time available to ensure the 
program is free of computational and logic errors, it cannot be 
considered validated. Any application of this program without 
additional verification is at the risk of the user. 
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A 

Area 

ft 2 
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acceleration 
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ft 
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sec 
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flight velocity 
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velocity 
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ft 
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I. INTRODUCTION 


A. GENERAL 

A FORTRAN computer program has been created to determine 
transient missile and aircraft mission aerodynamic heating 
parameters such as heat transfer coefficients, insulated skin 
temperatures and heat flow rates at some point on the surface 
a wing or a body of revolution. Required data supplied from 
the user include material properties and trajectory data such 
as Mach number, angle of attack, altitude, and time of flight. 
This program calculates aerodynamic heating parameters for 
subsonic, transonic, supersonic and hypersonic (no reacting 
chemistry) flight regimes. 

B. BASIC THEORY 

The governing equation which serves as the basis for the 
computer program incorporates both radiative and forced 
convective heat transfer processes. The thermal energy 
balance equation for this insulated skin heat transfer case is 
given by 

= hA(T aw - T e ) - oeA( T a * - T r 4 ) . (1) 

The equation written in this form assumes that the adiabatic 
wall temperature, T 1W is greater than the insulated skin 
temperature, T,, which is greater than the radiation reference 
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temperature, T r . The heat transfer coefficient for convection 
is obtained from the Nusselt number which is a function of the 
Reynolds and Prandtl numbers. These dimensionless numbers are 
defined as 

(a) Ngg = , (b) N pr = —, (c) N m = -jr • (2) 

Once the complex heat transfer coefficient calculation is 
performed and the other significant heat transfer parameters 
determined, a finite difference numerical method can be 
applied, yielding a skin temperature profile for the specified 
point of interest. 

C. SCOPE 

Chapter II details the theory used to develop the SKINTEMP 
program. 

Chapter III describes the actual computer program, its 
required entry data, expected output data and general program 
limitations. 

Chapter IV presents test case results with the appropriate 
comparisons, which include intermediate heat transfer 
calculations. One test case is obtained from a heat transfer 
textbook and compared to corresponding SKINTEMP results. A 
less descriptive NASA test case is also compared with SKINTEMP 
results. 

Chapter V provides recommendations for future code 
development. 


2 



II. TRANSIENT HEAT FLOW THEORY 


The fundamental concepts used to build the theoretical 
model for SKINTEMP involve two basic modes of heat transfer, 
convection and radiation. Forced convection is assumed because 
the vehicle is propelled through the air by the release of 
chemical energy. During supersonic and hypersonic flight, 
the local stagnation pressure at the edge of the boundary 
layer is assumed to be defined by the stagnation pressure 
behind a normal shock wave. Figure 1 illustrates this 
concept. 
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Boundary layer flow is assumed to be turbulent due to the 
magnitude of the Reynolds number which is generally greater 
than 500,000. The gas dynamic relations are based on inviscid 
flow, no reacting chemistry assumptions [Ref. 1]. 

The primary objective of the computer code is to rapidly 
compute the transient insulated skin temperatures along the 
trajectory of the aircraft at a given point somewhere on the 
surface. To accomplish this, the thermal energy balance 
equation 

mC p -|| = hA(T aw - T e ) ~ oeA(T e 1 - T*) (1) 

is solved for T, using a finite difference method [Ref. 2]. 
However, before this can be accomplished, the adiabatic wall 
temperature and the heat transfer coefficient must be 
determined. The required data for the skin temperature 
calculations include trajectory parameters (time, Mach 
number, angle of attack, altitude), material properties 
(density, material specific heat, emissivity), the location 
and skin thickness at the point of interest, and a radiation 
reference temperature. This radiation temperature is the 
temperature of a distant surface seen by the insulated skin 
element, which for this program, is either space or the earth. 
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A. THE ADIABATIC WALL TEMPERATURE 

The first parameter needed to determine skin temperature 
is the adiabatic wall temperature. The adiabatic wall 
temperature is given by 


Tjv = T l 


[l + 


(3) 


The local Mach number needed for the solution. 


Mr = 


P tL (T-D/r 


-1 


Y-l ' 


(4) 


is dependent on the local pressure and the local stagnation 
pressure. The speed of the aircraft, whether subsonic or 
supersonic will determine the local stagnation pressure. If 
the trajectory profile calls for subsonic flight, then the 
local stagnation pressure will be 


PtL - Pl 


1 + 


alT/T-i 


(5) 


If the Mach number is in the supersonic range, the local 
stagnation pressure is then 


PtL “ Pfr> 


2 


i + JL 

• 2 


r/r-i 


Y+l 


111 

Y + 1 


l/i-r 


where 


( 6 ) 


Pt- - 


- -P- [l + 




rt/r -1 


(7) 
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Freestream pressure, temperature, and density are calculated 
using standard atmosphere property value approximations found 
in Reference 3. 


Local surface pressure is determined using modified 
Newtonian theory, [Ref. 4] where 


p l = <Wcos 2 0g + p m 


~pmax 


P tl - P- 
Q 


q = l/2p.vf. 


( 8 ) 

(9) 

( 10 ) 


Modified Newtonian flow theory has been shown to be applicable 
for the prediction of local surface static pressures over all 
surfaces experiencing non-separated flow. [Ref. 5,6] The 
accuracy of Modified Newtonian flow theory is quite reliable 
at hypersonic speeds. Its use at subsonic and low supersonic 
speeds is justified by the continuity and simplicity it 
provides in the lower Mach number regions where aerodynamic 
heating rates are so low as to be negligible. 

The remaining parameter needed to determine the adiabatic 
wall temperature is the recovery factor (turbulent flow is 
assumed), given by 

i = Npg 2 (ID 


where the Prandtl number is 


(12) 

The specific heat of air used in the Prandtl number is given 
by [Ref. 2] 

C p = .219 + .342£-4T i - . 293E-87%. (13) 
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Thermal conductivity [Ref. 7] and dynamic viscosity [Ref. 2] 
are also temperature dependent and are given respectively by 


Jc = 


. 0141 


7 07 

Tl ' 

77+225 

L J 

492. 


2.21E-Q t\ 12 
7^+198.7 


B. TEE HEAT TRANSFER COEFFICIENT 


(14) 

(15) 


In order to calculate the heat transfer coefficient,one 
uses the concept of the reference temperature. The reference 
temperature is a temperature whose magnitude is typically 
between that of the wall temperature and freestream 
temperature and in low speed flow, the reference temperature 
does become the average of the two [Ref. 2]. This reference 
temperature more accurately models the temperature 
distribution that occurs within the boundary layer, which has 
a significant impact on the skin temperature [Ref. 8] . The 
transport properties should be evaluated at the reference 
temperature [Ref. 9] which is defined by 


T* = T l + .5 (T s - T l ) + .22 0V " • (16) 


Values for the thermal conductivity, dynamic viscosity, and 
specific heat used in the Nusselt number are based on this 
reference temperature. 

The Nusselt number can be defined as 

(17) 

k 
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Based on experimental data, the Nusselt number can also be 
written as a function of Reynolds and Prandtl number, [Ref. 2] 

Nm = .0292 (N^) s (Np R ) 1/3 (18) 

where the boundary layer is assumed to be turbulent and the 

Reynolds number is defined as 


Nre 


P 'Vi* 


(19) 


With these significant variables defined, the thermal 
energy balance equation can be solved using a simplified 
finite difference method. This method calculates the 
temperature value at a future time based on a calculated time 
increment and the current temperature value. 


C. THE SKIN TEMPERATURE 

Solution of the thermal energy balance equation for the 
incremental skin temperature yields 

r.. Jn = [ hA(T aw - T s ) - oeA ( T* - T r <) ] A* + T s<i . (20) 


For stability purposes, the time interval associated with this 
finite difference method is given by 


At £ 


mc p 

h+h z 


( 21 ) 


where 

h z = A€0[rj + t\T z + T 2 z T s + T z ] . (22) 


Each progressive skin temperature is dependent on the 
previous temperature, heat transfer coefficient and adiabatic 
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wall temperature. The time interval choice is critical and 
must meet the stability requirement of the finite difference 
method. This stability requirement, which prevents the 
solution from diverging, is given by 


At (h + H r ) 
mC p 


1 . 


(23) 


To meet this requirement, the time interval was chosen as 


At 


.25 


mC p 

h + H/ 


(24) 


D. THE HEAT FLOW RATE 


The heat flow rate is dependent on the location of the 
stagnation heating point. In the stagnation region for the 
axisymmetrical and two dimensional platforms, the maximum 
stagnation point is 70 degrees from the centerline [Ref. 10] 
and the stagnation heating rate is given by [Ref. 11] 


*stag 


_ 17600 


y/Rad N 


_e_ 

u 

3.15 

Hsl 

1 

03 

Po 

. u c. 


[”sl 

& *545J 


(25) 


Away from the stagnation point but prior to regions that can 
be treated as flat plates, the heat flow rate [Ref. 10] is 
assumed to be 


Qheat = QstagCOS (0) 3/2 . 


(26) 


Otherwise, in areas where flat plate assumptions hold, the 
heat flow rate is 

Qheat = hA(T aw - T s ) - oeA ( T s * - T t *) . (27) 
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III. THE SKINTEMP COMPUTER PROGRAM 


A. THE COMPUTER CODE STRUCTURE 

The flow logic used for the SKINTEMP program is 
illustrated in Figure 2. Th' program is divided into four 
major sections, data entry, atmospheric property generation, 
skin temperature calculations and data output. 

The first module contains the data entry loop. The data 
can be entered via two methods, keyboard entry or file entry. 
If the user chooses to utilize keyboard entry, the program 
will prompt the user for the specific data. Each entry should 
be followed by a keyboard hard return. The file entry method 
requires the user to create an external file using any 
standard editor containing all the required information. For 
file entry, the correct syntax to enter at the prompt is 
SKINTEMP.in. The .in refers to the entry file; any input file 
used must have the .in included in the file name. This 
command allows data to be read automatically by the input 
subroutine. 

The important note to remember is that for file entry, 
each piece of fixed point data must be placed on a separate 
line and the entry data must correlate with the data requested 
by the program. If the program terminates with an error 
message, this would be the first place to check for mistakes. 
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The second section of the code generates standard 
atmospheric values that are needed for the freestream 


pressure, temperature, and density calculations. The 
governing equations used to model atmospheric properties are 
divided into two regions, gradient and isothermal [Ref. 3]. 
The gradient relationships are 


p _ _T_-[g/LR* i] 
Po To 
p T ~ ^3/LR) 

~p~ ~ ~T 
c o o 

T = L(h-h 0 ) + T 0 . 


(28) 

(29) 

(30) 


There are several isothermal regions in the atmosphere. 
While the temperatures remain constant in each of these 
isothermal regions, the temperatures of the gradient regions 
will vary. [Ref. 3] The pressure and density relationships 
for isothermal regions are given by 


P _ -^h-h 0 ]/RT 

(31) 

Po 

_e_ = e -r h »] /iiT . 

Po 

(32) 


With these relationships, the atmosphere is modeled to 260,000 
feet with isothermal layers from approximately 36,000 to 
82,000 feet (temperature « 392.4 °R, density = .00071 slug/ft 3 , 
pressure = 476 lb/ft 2 ) and 154,000 to 174,000 feet 

(temperature = 509 °R, density = .00000303 slug/ft 3 , pressure 
= 2.6419 lb/ft 2 ). 






The third major section of SKINTEMP is the calculation of 
the skin temperature. This section is a large counted loop 
because each successive skin temperature is dependent on the 
previous temperature calculated. This loop is partitioned 
into four sections. The first section calculates adiabatic 
wall temperature and the second determines the heat transfer 
coefficient. These two calculations lead to the skin 
temperature solution, which is in the third section. The 
fourth and final section calculates the heat flow rate. [Ref. 
Chapter II] Appendix A contains the actual computer program. 
The initial documentation of the program contains a list of 
the variables, their definitions, and the English units used. 

It is important that the user understand the differences 
between the time values used throughout SKINTEMP. The user 
should input time intervals that best linearize the 
trajectory. In other words, whenever the Mach number, 
altitude or angle of attack is changing, this time and 
associated flight parameters should be a data entry point. 
These time intervals can be the print out time intervals but 
SKINTEMP also allows the user to print the results at 
different time intervals if desired. These input and output 
time intervals are independent of the time interval chosen for 
the heat balance calculations. Figure 3 demonstrates these 
differences and should be used as a reference when developing 
a data entry strategy for SKINTEMP. 
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BLOCK A 


READ INPUT DATA FROM USER OR FILE INTO ARRAY 


GENERATE ATMOSPHERIC TABLES 


BLOCK B 


TAW IS COMPUTED USING 

P t L' Qi 0£,/ Cpmax' ^L> ^L> 

T tL , k, C p , Ji, N pR , I, T l 

= I == 

h IS COMPUTED USING 

r r , P r ,|i r , v l , n &. 

Cp,k r ,Np R ,Nm 


id 


SKIN TEMPERATURE IS COMPUTED AFTER CALCULATING St 

t 

HEAT FLOW RATE COMPUTED BASED ON STAGNATION POINT 
- 1 - 

IF ^ 

£6 t < INPUT TIME INTERVAL 

_ GOTO BLOCK B _ 

I 

IF _ 

£5t ;> INPUT TIME INTERVAL 

PRINT OUTPUT 
GOTO BLOCK A 


Figure 2. SKINTEMP Flowchart 

See Chapter II for a complete description of the relationships 
needed to obtain the desired results. 
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TIME INTERVAL COMPARISONS 



INPUT TIME 


H—H— 1 -1 - I -4 - M -1—I—H—h- 

STEP TIME 


OUTPUT TIME 


Figure 3. Time Interval Comparison 
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B. REQUIRED DATA INPUT 

SKINTEMP allows the user to enter data through 
the keyboard or through an external file. In either case, the 
following data is required to successfully run the program: 

• The shape of the body, which must be either two 
dimensional or axisymmetric 

• If a two dimensional wing shape is selected, the wing 
sweep (degrees) 

• The launch surface skin temperature of the vehicle (°F) 


• For each time interval, the Mach number, the angle of 
attack (degrees), and the altitude (feet) 

*************************************** 

WARNING: Entering a Mach number of 
zero will cause the 
program to abort. See 
Section D. below. 

*************************************** 

• The distance from the nose or leading edge to the point of 
interest (feet) 

• The thickness of the surface at the point of interest 
(inches) 

• The nose radius (feet) 

• The material density (pound per cubic foot) 

• The material emissivity (dimensionless) 

• The material specific heat (BTU per pound-°F) 

• The orientation of the surface. Does the surface face the 
earth or space for the majority of the flight? This 
provides the radiation reference temperature of -460°F for 
space and 60 °F for the earth. (°F) 

• The angle from the unit surface normal of the skin element 
to the freestream (degrees) 

• The output time interval (seconds) 
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c. SKZxrrsMP output 


The output of the computer program is in tabular form. 
The user can choose any time interval for printing independent 
of the data time entries. See Appendix B for an sample 
output. SKINTEMP will also produce an external UNIX file 
named Fort.3 which produces a table of values for the time and 
skin temperature that can be used in typical engineering 
graphics/plotting programs. 

The program is clearly documented and the print subroutine 
is the last section in the program. The output format is as 
follows: 

TIME ADIABATIC WALL HEAT TRANSFER SKIN HEAT FLOW 

TEMPERATURE COEFFICIENT TEMPERATURE RATE 

(sec) (°Rankine) (BTU/sec-ft 2 -°F) (°Rankine) (BTU/sec) 

D. SKINTEMP LIMITATIONS 

Th3 purpose of this code is to provide a design tool for 
conceptual missile and aircraft design projects. At the 
present time, SKINTEMP can only calculate the transient heat 
parameters for bodies that are either two dimensional or 
axisymmetric body shaped. See Chapter V for future 
recommendations for code expansion. 

Another limitation and perhaps the most important for user 
information is that at no time can the Mach number entered be 
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equal to zero. The reason for this is that in the calculation 
for C r „... dynamic pressure is the divisor. C r „„ and dynamic 
pressure are given by 

c P ,max = ( p tL " PJ/Q (9) where q = i/2pv 2 = |pw 2 . (33) 

A Mach number of zero will cause C r „„ to be undefined. If 
this occurs, the program will terminate without results. If 
a Mach number of zero is the actual value achieved in the 
trajectory, a small value should be entered, i.e., .001 Mach. 

The last limitation is the array size of the input data. 
At the present time, the array size of data points entered is 
100. If more than 100 data points are desired by the user, 
the array size should be altered in the body of SKINTEMP. The 
array dimension statement is the first to occur after the 
initial documentation and the maximum limit is dictated by the 
version of FORTRAN being used. 




IV. TEST CASE RESULTS 


A. TRUITT TEST CASE 

The test case used to validate SKINTEMP is based on an 
example skin temperature problem from Robert W. Truitt's text 
Fundamentals of Aerodynamic Heating , [Ref 12} . The example 
missile profiled is axisymmetric in shape and 40 feet in 
length. The skin material used is Inconel steel and the skin 
thickness of the material at the point of interest, 20 feet 
from the nose, is .00504 inches. The flight trajectory data 
are given in Figure 4. 

Trajectory Mach number 




Figure 4. Truitt's trajectory data 













Truitt assumes turbulent flow but does not use specific 
local or reference state conditions. Instead, he assumes 
freestream conditions throughout the calculations. This 
approach is responsible for at least some of the noticeable 
deviation in the comparative results [Ref. Section 2]. 

A Runge - Kutta numerical integration method for the skin 
temperature calculation is employed by Truitt rather than the 
finite difference method used in SKINTEMP. This difference in 
approach was one of the primary reasons this test case was 
used. Similar results from two different equivalent numerical 
methods tend to confirm results. 

1. Data Comparisons 

Throughout Truitt's skin temperature calculation 
example several intermediate heat transfer and fluid flow 
parameter values are defined. These computed values include 
freestream static pressure and temperature, Reynolds number 
and the computed skin temperature (based on freestream 
properties). The comparisons of these parameters are listed 
in the following graphs and tables to validate SKINTEMP and 
illustrate some of the theoretical differences between 
calculation methods. 
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a. Freestream Pressure Comparison 

Figure 5 graphically displays the difference in 
freestream static pressures while Table 1 details the exact 
values computed. Truitt's values for freestream static 
atmospheric density are not as close to the NACA standard 
atmosphere [Ref. 13] as the SKINTEMP code calculations and the 
slight divergence in values contributes to the difference in 
the predicted skin temperatures. 


FREESTREAM PRESSURE COMPARISON 



5 10 IS 20 25 30 

Time (seconds) 


Figure 5. Freestream static pressure comparison 








TABLE 1. FREESTREAM PRESSURE COMPARISON FOR TRUITT’S 
TRAJECTORY 


TIME TRUITT SKINTEMP NACA 

(SECONDS) (LB/FT 3 ) (LB/FT 3 ) (LB/FT 3 ) 
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a. Temperature Comparison For Truitt's Trajectory 
Figure 6 and Table 2 compare Truitt's freestream 
temperature to SKINTEMP's freestream and reference static 
temperature. The freestream static temperature calculations 
compare favorably. Again, the divergence in values can be 
partially attributed to differences in atmospheric modelling. 

Local temperatures are compared with the 
freestream temperatures to illustrate the difference of using 
local atmospheric properties instead of freestream properties. 
SKINTEMP uses local conditions for the heat transfer 
calculations while Truitt utilizes the freestream conditions. 
Temperatures near the skin of a hypersonic vehicle will be 
much higher than the surrounding freestream air due to 
boundary layer interactions (i.e. skin friction). Thus, 
calculations using freestream properties are expected to be 
more approximate than those using the loaal conditions. 
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TABLE 2. TEMPERATURE COMPARISON FOR TRUITT'S TRAJECTORY. 



















































































































b. Reynolds Number Comparison For Truitt's 

Trajectory 

The difference in theoretical approaches is very 
apparent when comparing the Reynolds numbers in Figure 7 and 
Table 3. Truitt's Reynolds number calculations are very close 
to SKINTEMP's initial calculations, but again, the divergence 
in freestream atmospheric properties clearly has an impact on 
the freestream Reynolds number. The only significant 
consideration is that this difference becomes larger as the 
program builds on data from the previous time step. 

Figure 7 clearly illustrates that the reference 
state Reynolds number is significantly smaller than the 
freestream value. Recall that the reference state utilizes 
the local flow field values at the edge of the boundary layer 
with the transport properties evaluated at the reference 
temperature, and that the Reynolds number is given as 


n:, 


’RE 


p’VX 

n- ' 


where 


(18) 


p = P/RT 
V * MjyFF 
and 

2.21E-Q Tl /2 
T l + 198.7 


(34) 

(35) 

(15) 
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The Reynolds number is inversely proportional to the 
temperature. Use of the local values for temperature, which 
are significantly larger than the freestream values of 
temperature, will cause the Reynolds number to decrease. 

Another parameter that will contribute to a 
difference in Reynolds number calculations is the relationship 
used to determine viscosity, equation 15. Truitt did not 
utilize the same relationship and even though the values are 
of the same magnitude, the difference only adds to the 
deviation in results. 


FREESTREAM AND REFERENCE REYNOLDS NUMBER COMPARISON 



Figure 7. Reynolds Number Comparison 
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TABLE 3. REYNOLDS NUMBER COMPARISON FOR TRUITT'S TRAJECTORY 
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Skin Temperature (Rankine) 


c. Skin Temperature Comparison 

Figure 8 and Table 4 compare the final skin 
temperature calculations between Truitt and SKINTEMP. 


2500 


TRUITT SKIN TEMPERATURE PREDICTION COMPARISON 
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Figure 8. Skin Temperature Comparison 








TABLE 4. SKIN TEMPERATURE COMPARISON FOR TRUITT'S TRAJECTORY 


TIME 

(SECONDS) 


TRUITT 

(°RANKINE) 


1740 




2150 




2250 


2250 


2250 


2250 


SKINTEMP 

(°RANKINE) 



1091 


14 


1761 


1917 


1845 


1781 


1711 


1626 


1704 


1968 


2132 


225 


2375 


231 


2297 


2300 


2310 


2071 
























































































1. Results Prom Truitt's Trajectory 


The graphs and tables Indicate comparable results in 
the two different methods of calculating skin temperature. 
There are numerous reasons for the slight differences in the 
results. The discrepancy in standard atmosphere values leads 
to significant differences in the aerodynamic parameters, 
which cause marked deviations in the heat transfer 
coefficients, and the adiabatic wall and skin temperatures. 

Another major difference is the use of the reference 
state in SKINTEMP and the freestream state in Truitt's test 
case. Using the reference condition instead of the freestream 
state is considered to be the more correct model of what is 
actually occurring aerodynamically. One aspect of this 
difference is manifested in the Reynolds number calculation. 
Despite the significant difference in Reynolds numbers used to 
compute the Nusselt number, the results are nearly the same. 
Recall that to compute the Nusselt number, the Reynolds number 
is raised to the 0.8 power. Large differences in Reynolds 
number between the two methods are minimized in this 
calculation, allowing the skin temperatures to more closely 
correspond with one another. 

Truitt also uses varying skin material specific heats 
while SKINTEMP assumes an average constant specific heat. 
This does not cause a noticeable change in skin temperature 
calculations but does account for a few degrees of difference. 
The choice of time constants is also critical for stability 




purposes. Truitt uses a time step of two seconds while the 
average time constant step generated by SKINTEMP is 0.3 
seconds. As indicated by Truitt, the smaller the time 
constant, the more accurate the results. Truitt used the 
larger time constant for ease in calculations because they 
were all hand computed. With all other variables being equal, 
the calculation of skin temperatures with the more simplified 
finite difference method versus the Runge-Kutta method is 
expected to give the same degree of accuracy, provided the 
time step is small. 

While the Runge-Kutta integration method is an 
accurate numerical method, Truitt did not have the same 
computer capabilities that are in use today and the current 
standard atmosphere tables are considered to be more accurate 
than those Truitt used. Even though the SKINTEMP results are 
on average nine percent greater than Truittt's results, the 
skin temperature curves are similar with SKINTEMP producing 
the more conservative worst case results. 

A. NASA TEST CASE 

NASA AMES provided an additional test case for SKINTEMP 
verification. The trajectory is typical of a Single-Stage - 
To-Orbit (SSTO) rocket powered vehicle and the point of 
interest is 75 feet aft of the nose. The thickness of this 
point of interest is 0.035 inches and the material at this 
point is Tailorable Advanced Blanket Insulation (TABI). 
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There is a body angle of incidence of 13 degrees in addition 
to the varying angle of attack experienced through flight. 
Figure 9 graphically depicts the results and is followed by 
the detailed listing of skin temperatures. 


NASA SKIN TEMPERATURE PREDICTION COMPARISON 



Figure 9. NASA Skin Temperature Comparison 
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TABLE 5. SKIN TEMPERATURE COMPARISON FOR NASA AMES TEST CASE 


TIME 

(SECONDS) 


2.61 


3.48 


4.69 


5.78 


6.82 


7.70 


8.71 


9.81 


10.76 


11.79 


12.85 


13.72 


14.72 


15.75 


16.75 


17.75 


18.75 


19.70 


20.78 


21.38 


NASA 

(°RANKINE) 

SKINTEMP 

(°RANKINE) 

511 

511 

571 

570 

590 

630 

738 

870 

1313 

1104 

1570 

1449 

1821 

1756 

2094 

1979 

2318 

2262 

2463 

2424 

2637 

2550 

2799 

2645 

2912 

2733 

3011 

2806 

3127 

2864 

3200 

2911 

3265 

2974 

3308 

3019 

3267 

3084 

3008 

2606 

3231 

2133 
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1. Results From NASA Trajectory 

Again, the results from the two skin temperature 
calculation methods compare favorably. The NASA test case was 
also based on freestream conditions with surface pressures 
given by tangent wedge/cone theory. While the use of 
freestream conditions is not theoretically correct, there has 
been significant experimental data generated to justify this 
assumption [Ref. 14] . However, these assumptions are best 
used for upper surfaces. SKINTEMP is not limited to this 
assumption and can correctly calculate skin temperatures for 
upper and lower surfaces. 
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V. CONCLUSIONS AND RECOMMENDATIONS 


A. CONCLUSIONS 

The SKINTEMP computer program was developed as a tool for 
use in the conceptual design of high speed missiles and 
aircraft. It can be used over a wide range of trajectories, 
profiles and configurations and rapidly yields results that 
enable the user to easily investigate changes in material 
characteristics and flight parameters without having to employ 
more complex and cumbersome calculations. As with any 
computer program, the results can only be as good as the data 
entered. This is especially true of SKINTEMP because valid 
results can only be expected if the trajectory is accurately 
modeled. 

Validation of SKINTEMP was achieved using a textbook skin 
temperature calculation example and a more current NASA test 
case. SKINTEMP compared with the NASA and Truitt test cases 
suggest that the SKINTEMP results are reasonable. While each 
of the methods has its own merits, the fundamental theoretical 
base, direct approach and simplicity of SKINTEMP make it ideal 
as an aerodynamic and astrodynamic design tool. 








B. RECOMMENDATIONS 


SKINTEMP currently calculates the insulated skin 
temperature and heat flow rate at only one surface location. 
The computer code could be expanded to perform a thermal nodal 
analysis that could predict skin temperatures along an 
extended surface. Moreover, temperatures of materials within 
the air/spacecraft affected by the skin temperature could be 
calculated. This would be very useful in the design process 
to initially identify hot spots which, especially in 
hypersonic flight, could be crucial to the overall design. 
SKINTEMP was written strictly for the turbulent flow case 
because this generally produces worst case skin temperatures. 
The code could easily be modified to model laminar flow if 
turbulent flow was not the required assumption. 

Because the skin temperature problem is quite detailed and 
the inputs are extensive, complete trajectory and 
corresponding surface skin temperatures profiles are difficult 
to find in the literature. Additional test cases could be 
developed for SKINTEMP for further validation as well as any 
experimental data that is available. 

The important objective to remember in any significant 
code modification of SKINTEMP is that the program needs to run 
quickly and allow the user to modify any design parameters. 
If the program becomes too complex and time consuming, then 
the usefulness of SKINTEMP is degraded. 
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APPENDIX A. SKINTSMP COMPUTER CODE 


LT Nary L* Cummings 
591 

This program calculataa transient aisilon aerodynamic heating 
parameters (Adiabatic vail temperatures, local heat transfer 
coeffcients, and insulated skin temperatures) using various 
basic flight parameters as input. 


Variable declarations 

BODY - determines shape of body, either wing or 
axisymmetric 

SWBEP • angle of wing sweep (degrees) 

TEMP - launch vehicle temperature (Fahrenheit) 

NUMDAT - number of lines of data 

TABLE • input data consisting of time (sec), Mach number, 
of attack (degrees), and height (feet) 

X - distance from leading edge to point of Interest (feet) 

THIN - thickness of interest point (Inches) 

RAD - nose radius of body (ft) 

MATDENS - material density (lbm/ft'3) 

MATEMIS - material emissivity 

MATSPHT - material specific heat <BTU/lbm-r) 

TEMPRAD - radiation temperature (rahrennheit) 

ANGLE • angle from the unit surface normal to the 
freestream (degrees) 

PRINTOUT - user defined print time (seconds) 
g - acceleration/weight constant (Ibm-ft/lbf-eec~2) 

LI - temperature lapse rate from 0 - 36009 ft (R/ft) 

L2 * temperature lapse rate from 62021 - 154200 ft (R/ft) 

R - gas constant (ft-lbf/slug-R) 

TH - thickness of interest point (feet) 

Densin - initial density at sea level (slug/ ft~3) 

Pin - initial pressure at sea level (lbf/ft*2) 

Tempin - launch vehicle temperature (Rankine) 

T/atmT - freestream temperature (Rankine) 

P/atmpr - freestream pressure (lbf/ft~2) 
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DENS/atmdena - freeatream danaity (alug/ ft*3) 

PreaaT - praaaura at tha tropoephere calling (lbf/ft~2) 
DanaT - danaity at tha tropoephere calling (alug/ ft*3) 
Ptl - local atagnation praaaura (lbf/ft*2> 

Pt - freeatream atagnation praaaura (lbf/ft‘2) 

Tt * freeatream atagnation taaiparatura (Rankine) 

Ma - calculatad fraaatraan Mach number 
Ml • Hach number behind ehock wave 
Hi - local Hach number 
a - fraaatream aonic velocity (ft/aec) 

V - fraaatream velocity (ft/aec) 
q - dynamic praaaura (alug/ft-aec"2) 

gamma - ratio of apeciflc heat at conatant praaaura to 
apacific heat at conatant volume 

aigma - Stefan-Boltzmann conatant (BTU/hr-ft*2-R A 4) 

theta - angle of attack added to the normal fraaatream 
angle (radiana) 

CpMax - Maximum praaaura coefficient 
Pi - local praaaura (lbf/ft~2) 

Tt1 - local atagnation temperature (Rankine) 

Tl - local temperature (Rankine) 

Tn - akin tampertaura (Rankine) 

TAW - adiabatic wall temperature (Rankine) 

re - recovery factor 

Pr - Prandtl number 

u • viacoaity (alug/aec-ft) 

k - thermal conductivity (BTU/aec-ft-r) 

Cl - Sutherland'a conatant (slug/ft-aec-R~.5) 

C2 - Sutherland'a conatant (Rankine) 

Cp « apacific heat (BTU/lbm-R) 

Tatar - reference temperature (Rankine) 
denatar - reference danaity <alug/ft*3) 
uatar - reference viacoaity (alug/aec-ft) 
al - local aonic velocity (ft/aec) 
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VI - local velocity (ft/sec) 

Restar - reference Raynolda number 

Cpatar - reference specific hast (BTU/lba-R) 

kstar - rafaranca thermal conductivity (BTU/sec-ft-P) 

Prstar - rafaranca Prandtl number 

NUstar - reference Nusselt number 

h * local heat| transfer coefficient (BTU/sec-ft*2-F) 
w - specific weight (lbf/ft*3) 

TEHPNEX - skin temperature calculated based on the 
previous skin temperature (Rankine) 

Hr - tiae differential calculation aid 

tiaa - time step (sec) 

suatiaa * interval tiae counter (sec) 

runtime • total tiae counter (sec) 

printtiaa - tiae printing counter (sec) 

stagpt - stagnation point (ft) 

ANGLER - Angle of attack (radians) 

hsi - stagnation enthalpy (BTU/lba) 

hws - wall enthalpy (BTU/lba) 

hw540 - wall enthalpy at 540 U (BTU/lba) 

qstag • stagnation point heating rate (BTU/sec) 

qheat - heat transferred (BTU/sec) 


dimension TABLE(100,4), T(100), p(100), DENS(IOO) 
integer BODY , NUHDAT 

real SWEEP, TEMP, TABLE, X, THIN, RAD, MATDENS, MATEMIS, 
C MATSPHT, TEMPRAD, ANGLE, TH, PRINTOUT, g, Ll, L2, L3, R, 
c Densln, Pin, Tenpin, T, P, DENS , PressT, DeriBT, 
c atapr, atadens, atmT, Pti, Pt, Tt, Ml, a, V, q, Ma, 
c gaaaa, sigma, theta, pi, printtiae, 

c CpMax, Pi, Ml, Ttl, Tl, Tn, TAW, re, Pr, u, k. Cl, C2, 
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c Cp, Tatar, denatar, uatar, al, VI. Mttac. Cpstar, 
c katar, Pcatar, NUatar, b. w, TEMPNEX, hr, tine, 
c atagpt, ANGLER, hal, hwa, hw540, qatag, qhaat, 
c tiaedif, auatiaa, daltatiaa, daltaa 


Input 


print * 

print *, 'If tha body ia axiaynaetric, antar 1.' 
print *, 'If tha body ia a wing, antar 2.' 

raad *, BODY 

print * 

if (BODY .EQ. 2) than 

print *, 'What ia tha wing aweap?' 
raad *, SWEEP 
print * 
and if 


print *, 'What ia tha launch taaparatura of tha vehicle? (r)‘ 
read *, TEMP 
print * 

print *, 'How aany data pointa will be entered for tine, nach nun 
c bar, angle of attack,' 
print *, 'and altitude?' 

raad *, NUMDAT 

print * 

print *, 'Enter tha tine (a), Mach nuaber, angle of attack (degr 
c aaa), and altitude (ft) in a row/colunn format.' 
print *, 'Hit return after every entry.' 
print * 

print *, '***WARNING*** Do not uae any data linea with H - 0' 

Do 20 ROW-1, NUHDAT 
Do 10 COL-1,4 

read *, TABLE(ROW,COL) 

10 continue 

20 continue 
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print * 

print * f 'What is the distance from the tip of the vehicle to the 
c point of interest? (ft) 1 

read *» X 

print * 

print *, 'What is the thickness of the surface at the point of in 
c terest? (in)' 

read *, THIN 

print * 

print *, 'What is the nose radius? (ft)' 
read *, RAD 
print * 

print *, 'What is the density of the Material? (lbm/ft*l)' 
read *, MATDENS 
print * 

print *, 'What is the emissivity of the Material?' 
read *, MATEMIS 
print * i 

print *, 'What is the specific nt-_ jf the material used? (BTU/lb 
cM-r)' 

read *, MATSPHT 
print * 

print *, 'For the majority of the flight path, if the surface rad 
elates towards the earth, enter 60.' 

print *, 'If it radiates to space, enter 0.' 

read *, temprad 

print * 

print *, 'What is the angle from the unit surface normal to the 
cfreestream? (degrees)' 

read *, ANGLE 

print * 

print *, 'What time increment would you like the skin temperature 
cs printed? (seconds)' 

read *, PRINTOUT 

print * 
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so 


gaaaa - 1.4 

slgaa - .1714E-8/3600.0 
pi - 3.1415927 
TH - THIN/12.0 
suatiaa - 0. 
daltatiaa - 0. 
printtiae - printout 

It (TEMPRAD .EQ. 60) then 

TEHPRAO - TEMPRAD + 460. 

end if 


do 100 n - 1, NUNOAT 

if (n .EQ. 1) then 
Tn - Teapin 
Ma - TABLE!n,2) 
atapr - P(n) 
atadens - Dens(n) 
ataT - T(n) 

if (BODY .EQ. 2) then 

Ma ■ Up * cos(sweep * pi/180.) 
end if 

tiaedif - (TABLE!n+1,1) - TABLE!n>1)) 
GOTO 60 
end if 

deltaa - (TABLE!n+1,2) - TABLE!n,2)) 
tiaedif - ( TABLE(n+1,1 ) - table( n,1)) 

if (euatiae .GT. 0.) then 
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* 

* 


Ataoapharic data calculation* 


* 

* 


9 - 32.174 
LI - -.0035665 
L2 - .001645 
L3 - -.002469 
R - 1716. 

Densln - .002378 
Pin - 2116.0 
Tenpin - TEMP ♦ 460.0 

do 30 HGT - 1, NUHDAT 

If (TABLE(HGT,4) .LE. 36089.0) then 
Troposphere 

T(HGT) - Tenpin + Ll * TABLE(HGT,4) 

P(HGT) - Pin * (T(HGT)/Tenpln)**(-g/(Ll * R )) 

DENS(HGT) - Oensln * ((T(HGT)/Tenpln)**-(g/(Ll * R) 

c ♦ m 

goto 30 
end If 

Stratosphere 

If (TABLE(HGT,4) .GE. 82021) then 

If (TABLE(HGT,4) .GE. 173885) then 

T(HGT) - 509 + L3 *(TABLE(HGT,4) - 173885) 

P(HGT) - 1.5 * (T(HGT)/509)*M-g/(L3 * R )) 

DENS(HGT) - .0000014 * ((T(HGT)/509)** 

C -((9/(L3 * R))+ 1.0)) 

goto 30 

end If 

If (TABLE(HGT,4) .GE. 154200) then 
T(HGT) - 509 

PressT - 2.6419 / 
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DanaT - .00000303 

P(HGT) - PraaaT * EXP (-(ig / U * T(HGT))) * 
C ( TABLE(HGT, 4) - 154200.))) 

DENS(HGT) - DanaT * P(HGT) / PraaaT 

goto 30 

•nd if 

T(HGT) - 389.9 t U * (TABLE(HGT,4) - 82021.) 

P(HGT) - 55.8 * (T(HGT)/389.97)**(-g/(L2 * R )) 

DENS (HGT) - .0000789 * ((T(HGT)/389.97)* * 

C -<(g/(L2 * E) )+ 1.0/) 

goto 30 

end if 

T(HGT) - -67 < + 460.0 
PreasT - 476 
DanaT - .000710 

p(HGT) - PraaaT * exp (-((g / (a * t(hgt))) * ( 

C TABLE ( HGT, 4) - 36089.))) 

DENS(HGT) - DanaT * P(HGT) / PraaaT 
30 continue 


Header 


print * 
print * 
print * 
print * 


print 40, 'Tina', 'Adiabatic Hall', 'Heat Tranafer', 
c'Skin Teap', 'Heat Tranafer Hate ' 


print 45, '(aec)', ' Teaperature (a)', 
c ' (BTU/aec-one aq ft)' 

print a 


' Coefficient', 
I 


(a)', 


40 

roraat(2X, 

A4, 

2X, 

A16, 

IX, 

A16, 3X, 

A9, 

5X, A21) 

45 

Poraat(2X, 

A5, 

2X, 

A16, 

4X, 

All, 6X, 

A4, 

8X, A21) 



* Subaonio va auperaonic calculationa * 

* * 
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H» - auktiae * deltea/(tiaedif} + TABLE(n,2) 


atape - P(n) + suatiae * (P(n+1) - p{n))/tiaedi£ 

atadens - Dens(n) + suatiae * (Dens(n+1) - Dens(n)) 
c /tiaedif 

ataT - T(n) + suatiae*(T(n+l)-T<n))/tiaedif 

if (BODY .EQ. 2) then 

Ha - Ha * cos(sweep * pi/180.) 
end if 

end if 

if (suatiae .EQ. 0 ) then 

Ha * TABLE( n , 2 ) 
atapr - P(n) 
atadens - Deas(n) 
ataT - T(n) 
if (BODY .EQ. 2) then 

Ha-Ha * cos(sweep * pi/180.) 
end if 

end if 


60 if (Ha .LE. 1.0) then 

Ptl - atapr * (1.0 + ((gamma-1.0)/2.0) * Ma**2)** 
c (gamma/(gamma-1.0)) 

else 

Pt - atapr * (1.0 + ((gaaaa-1.0)/2.0) * Ha**2)** 
c (gamma/(gamma-1.0)) 

Tt - ataT * (1.0 + ((gamma-1.0)/2.0) * Ha**2) 

Ptl - Pt * (( (gaama+1.0)/2 * Ma**2)/(1 + (gamma 
c ~1.0)/2 * Ma**2))**(gamma/(garoma-1.0)) * 

c (2.0*gaama/(gaaaa+1.0) * Ma**2 -(gamma-1.0)/ 

c (gamma+1.0))**(l/(1-gamma)) 


end if 


* * 

* Adiabatic wall temperature calculations * 

* * 





a - (gamaa * R * ataT)**.5 
V - a * Ha 

q - .5 * atadens * V**2 

theta - (ANGLE + TABLE(n« 3)) * pi/180. 

CpHax - (Ptl - atapr)/q 

Pi - q*CpHax * (costtheta))**2 + atapr 

Hi - (((Ptl/Pl)**((gaaaa-1.0)/gaaaa)-l.0) * (2.0/ 
c (gaaaa-1.0)))**.5 

Ttl - ataT * (1.0 + (gaaaa-1.0)/2.0 * Ma* A 2) 

Tl - Ttl/( 1.0 + (gaaaa-1.0)/2.0 * Ml**2) 

k - (.0141 * (717.0/(Tl + 225.)) * (Tl/492.0) 
c **1.5)/3600.0 

Cp - .219 + .342E-4 * Tl - .293E-8 * Tl**2 

Cl - 2.27E-8 

C2 - 198.7 

u - Cl * Tl**l.5/(Tl + C2) 

Pr - Cp * u * g /k 
ce - Pc**(l.0/3.0) 

TAW - Tl * (1.0+re * (gaama-1.0)/2.0 * Ml**2) 


Heat transfer coefficient calculations 


Tstar - Tl + .5 * (Tn-Tl) + .22 * (TAW-Tl) 

denstar * Pl/(R * Tstar) 

ustar - (Cl * Tstar**1.5)/(Tstar + C2) 

al - (gaaaa * R * Tl)**.5 

VI - al * Ml 

Restar - denstar * vl * X/ustar 

Cpstar - .219 + .342E-4 * Tstar - .293E-8 * Tstar 
c * *2 

kstar - (.0141 * (717.0/(Tstar + 225.)) * (Tstar/ 
c 492.0)**1.5)/3600.0 

Pretar - Cpstar * g * ustar/kstar 
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c 


c 

c 


if (Raster .LS. 1.E7) then 

NUstar - .0292 * (Raster**.8) * (Prstar 

**( 1 . 0 / 3 . 0 )) 

• Isa 

NUstar - .185 * (Restar) * (Prstar)** 

(1.0/3.0)/((loglO(Restar))** 
2.584) 

end if 

h • NUstar * kstar/x 



w - MATDENS / 9 

Hr - HATEHIS * sigma * (Tn**3 + Tn**2 * TEMP 
C RAO + TEMPRAD**2 * Tn + TEHPRAD* * 3 ) 

time - .25 * MATSPHT * w * 9 * th /(hr + h) 

TENPNEX - time*(((h * 1 *(TAW-Tn)) 
c - sigma * MATEMIS * (Tn»*4-Tr**4))/ 

C (MATSPHT *W*g*TH*l))+Tn 


* Stagnation point heat transfer calculation 


ANGLER - ANGLE * pi/180. 

hsl - Cp * Ttl 

hws - Cp * Tn 

hw540 - Cp * 540. 

stagpt - cos(70 * (pi/180.)) 

if (ANGLER .LE. stagpt) then 


qstag - 17600.0/SQRT(RAD) * SQRT(atmden6/Densin) * (V/ 
c 26000.0)**3.15 * ((hsl-hws)/(hsl-hw540)) 


if (BODY .EQ. 1) then 

qheat - qstag * (cos(ANGLER)**(3.0/2.0)) 
else 

qheat - qstag * (cos(ANGLER)**(3.0/2.C))/SQRT(2.0) 
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end if 


Heat transfer calculation 


qhaat - h * 1 * (TAW-Tn) - eigaa * 

C If ATEN IS * (Tn**4-TEHPHAD**4) 

and if 

Tn - Teapnex 

auatiaa - suatima + tiaa 
runtiaa - runtime + tiaa 
Tf - Taapnax 

if (auatiaa .GE. tiaedif) than 
daltatiae - auatiaa - tiaedif 

if (runtiaa .GE. printtiaa) then 

Print 80, printtiaa, TAW, h, Tf, qhaat 
WRITE(3,90) printtiaa, TEMPNEX 
printtiaa - (printtiaa + printout) 
and if 
auatiaa - 0 
goto 100 
else 

if (n .GE. nuadat) then 
goto 100 
end if 
GOTO 50 
end if 


* * 

* Output * 

* * 
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roraatt ox, n.2 t sx, r9.2, lox, re.6, 5x, r9.2, I2x f re.3) 


80 

90 Poraat(IX, r7.2, 10X, r9.3) 


100 continue 

print * 
print * 
print * 


end 


I 
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APPENDIX B. EXAMPLE SKINTEMP INPUT AND OUTPUT 


What l«t the launch temperature pf l lie vehicle? (F) 

Now m#n\ d'it« point! will he entered for time. Hitch aum her, angle of attach, 
and altitude? 

Enter the tint* fa), H»cH number, Angle oC attack (dogr ecsl # and altitude (ft l 
In a row/colwait format. 

Nit return after every entry. 

* 4 - WNNM?MO* • * Do not use any data line* with M « 0 

What h the distance from the tip of the vehicle to the point of Interest? (ft) 
What la the thickness of the surface at the point of in terest? (In) 

What t* the nose ladtus? <(t.l ‘ 

What Is the density of the material? (1bm/ft*3) 

What I* the emlsslvlty of the material? 

What Is the specific heat of the material used? (BTU/Ihm-D 

for the majority of the flight pal h, If the surface radiates towards the earth, 
enter *0. 

Ik’ 11 rad! ales to space, enter 0. 

What Is the angle from the unlf surface normal to the freest ream? (degrees) 
What time Increment would you like the skin temperattires printed? (seconds) 

SKINTEMP data entry prompts ™~ 
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1 

51 

32 

1 

.34 

0 

2*0 

2 

.718 

0 

500 

3 

1.08 

0 

750 

4 

1.445 

0 

1000 

5 

1.814 

0 

1500 

6 

1.735 

0 

2000 

7 

1.45 

0 

3000 

• 

1.570 

0 

5000 

9 

1.48 

0 

9500 

10 

1.397 

0 

10000 

11 

2.107 

0 

13000 

12 

2.83 

0 

15000 

13 

3.575 

0 

18000 

14 

4.348 

0 

22000 

15 

5.157 

0 

25000 

14 

5.095 

0 

27000 

17 

5.03 

0 


Date entered using the file entry method (Truitt test case) 
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34000 

It 

4.334 

0 

37000 

13 

4.104 

0 

40000 

20 

4.443 

0 

44000 

21 

3.30 

0 

30000 

22 

4.331 

0 

52000 

23 

i. 103 
0 

42000 

24 

3.254 

0 

44000 

25 

10.373 

0 

43000 

24 

11.521 

0 

•5000 

27 

12.47 

0 

•3500 

23 

13.317 

0 

35000 

23 

14.43 

0 

117000 

30 
15 
0 

123000 

31 
15 
0 

125000 

32 
15 
0 

125000 

20 

.005004 

.05 

530 

.3 

.121 

0 

30 

1 

1 


Data entered using the file entry method (Truitt test case, continued) 


53 







Tim*! 

Adiabatic Wall 

Heat Transfer 

Skin Temp 

Heat Transfer Rate 

(sec) 

Tomfifsrature (0) 

Coefficient 

(R> 

(BTU/sec-one sq ft) 

1.00 

529.30 

0.008575 

522.02 

0.041 

2.00 

609.95 

0.019815 

560.91 

1.244 

3.00 

698.44 

0.024123 

633.11 

2.018 

1.00 

793.95 

0.025823 

717.28 

2.503 

5.00 

794.17 

0.024537 

772.03 

0.526 

6.00 

760.86 

0.023045 

767.06 

-0.392 

7.00 

725.14 

0.02118G 

742.17 

-0.660 

? .00 

674.31 

0.017698 

704.24 

-0.854 

7,00 

656.61 

0.017273 

677.72 

-0.611 

»0.00 

792,45 

0.019156 

707.28 

2.048 

11.00 

1012.91 

0.019434 

850.20 

4.754 

12.00 

1397.34 

0.017993 

1091.53 

6.712 

13.00 

1859.94 

0.015438 

1408.19 

7.561 

11.00 

2451.81 

0.013316 

1761.73 

7.877 

15.00 

2457.28 

0.011G19 

1917.23 

0.815 

16.00 

2257.59 

0.006788 

1845.07 

-2.340 

17.00 

2202.84 

0.008247 

1781.19 

-1.444 

18.00 

2096.57 

0.007287 

1711.37 

-1.503 

13.00 

1979.69 

0.005898 

1626.52 

-1.607 

70.00 

2740.36 

0.005182 

1704.54 

3.274 

71.00 

3870.34 

0.004645 

1967.57 

5.554 

77.00 

5010.18 

0.003351 

2131.51 

2.475 

23.00 

6573.34 

0.002756 

2259.56 

2.906 

24.00 

8142.86 

0.002458 

2374.87 

2.701 

25.00 

9968.20 

0.001562 

2319.48 

-2.354 

26.00 

12715.81 

0 001156 

2297.95 

0.757 

o 

o 

12728.18 

0.001152 

2300.75 

0.071 

28.00 

15484.49 

0.000928 

2310.05 

0.229 

33.00 

18435.02 

0.000420 

2071.97 

-5.430 

70.00 

20525.59 

0.000321 

1952.78 

-1.970 

71 .00 

20525.59 

0.000323 

1936.63 

-0.227 

77.00 

70525.59 

1 

0.000324 

1935.69 

-0.013 


SKINTEMP output (Truitt test case) 
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